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Abstract 



Transition state theory was introduced in the 1930s to account for chemical reactions. Central to this theory is the idea of a potential 
{N| energy surface (PES). It was assumed that such a surface could be constructed using eigensolutions of the Schrodinger equation for 
,— I the molecular (Coulomb) Hamiltonian but at that time such calculations were not posssible. Nowadays quantum mechanical ab- 



initio electronic structure calculations are routine and from their results PESs can be constructed which are believed to approximate 



those assumed derivable from the eigensolutions. It is argued here that this belief is unfounded. It is suggested that the potential 
energy surface construction is more appropriately regarded as a legitimate and effective modification of quantum mechanics for 



C/3 



chemical purposes. 



in ■ 

; 1. Introduction 

<— | The principal aim of much of contemporary quantum chem- 
Q_|ical calculation is to calculate a potential energy surface from 
solutions of the Schrodinger equation for the clamped-nuclei 
Ch electronic Hamiltonian which gives the electronic energy at a 



choice of nuclear positions. These electronic energies are then 



qH added to the clamped-nuclei classical Coulomb repulsion en- 
— 1 ergy and the resulting "total" energies fitted to a functional form 
to provide a surface which, for a system of A > 2 nuclei, is of 
dimension 3A - 6. Nuclear motion is treated as occurring on 
this surface. 

The idea of a potential energy surface in such a role began 



> 

On 



with the work of Eyring and the almost contemporary work of 
^J. Polanyi in the 1930s. Its basis is described in Chapter 16 of the 
\Q textbook by Eyring, Walter and Kimball, published in 1944 |Q]]. 



O 
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When describing the reaction between the hydrogen molecule 
and a hydrogen atom, the authors say: 

The properties of this system are completely de- 
termined by the Schrodinger equation governing it, 
that is by its eigenfunctions. The system itself may 
be represented by a point in four dimensional space: 
three dimensions are required to express the relative 
positions of the nuclei, one additional dimension is 
required to specify the energy. (We assume here that 
the motion of the electrons is so rapid that the elec- 
trons form a static field for the slower nuclear mo- 
tions.) The three internuclear distances... can con- 
veniently be used to specify the configuration of the 
system. 

They then give the name "potential energy surface" to the pre- 
sumed calculated entity and go on to explain how transition 
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probabilities for the hydrogen exchange reaction can be calcu- 
lated using its characteristic shape. They take it as obvious that 
the nuclei can be treated as classical particles whose positions 
can be fixed and do not mention the work of Born and Oppen- 
heimer which, 17 years before the publication of their book, 
had argued that such a treatment of the nuclei could be justified 
if the potential energy surface had a unique minimum and that 
the nuclear motion involved only small departures from equi- 
librium |H H . 

In the textbook by Pauling and Wilson J3l the work of Born 
and Oppenheimer is quoted as justifying their use of a wave- 
function which is a simple product of an electronic and a nu- 
clear part in order to describe the vibration and rotation of 
molecules. They then introduce the potential energy function, 
which they call the "electronic energy function", for a diatomic 
system. In their description of a chemical reaction involving 
three atoms they attribute the idea of fixing the nuclei to the 
work of London who used such an approach in his 1928 pa- 
per fl. It is interesting to note that London actually called 
this approach "adiabatic" saying that he assumed that as the nu- 
clei moved they acted as adiabatic parameters in the electronic 
wavefunction. 

Although treating the nuclei as classical particles that may be 
fixed in space to generate a potential energy surface is often now 
referred to as "making the Born-Oppenheimer approximation" 
it is something of a mis-attribution as can be seen from this 
discussion. This was indeed recognized by Born and in 1951 
he published a paper (6J in which he presented a more general 
account of the separation of electronic and nuclear motion than 
that originally offered. It is to this later work by Born that the 
more general idea of a potential energy surface is regarded as 
owing its justification and it is that work that will be considered 
here. 

We aim to show that the potential energy surface does not 
arise naturally from the solution of the Schrodinger equation 
for the molecular Coulomb Hamiltonian; rather its appearance 
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requires the additional assumption that the nuclei can at first be 
treated as classical distinguishable particles and only later (af- 
ter the potential energy surface has materialized) as quantum 
particles. In our view this assumption, although often very suc- 
cessful in practice, is ad hoc. 

The outline of the paper is as follows. In the next section 
we review the main features of the standard Born-Oppenheimer 
and Born adiabatic treatments following Born and Huang's 
well-known book |7|]. The key idea is that the nuclear kinetic 
energy contribution can be treated as a small perturbation of the 
electronic energy; the small parameter (k) in the formalism is 
obtained from the ratio of the electronic mass (m) to the nuclear 
mass (M ): m/M = k 4 . The argument leads to the expres- 
sion of the molecular Hamiltonian as the sum of the "clamped- 
nuclei" electronic Hamiltonian (independent of k) and the nu- 
clear kinetic energy operator (oc k 4 ), equation ©■ 

In $3]we attempt a careful reformulation of the conventional 
Born-Oppenheimer argument drawing on results from the mod- 
ern mathematical literature. The calculation is essentially con- 
cerned with the internal motion of the electrons and nuclei so 
we require the part of the molecular Hamiltonian that remains 
after the center-of-mass contribution has been removed. We 
show that it is possible to express the internal motion Hamilto- 
nian in a form analogous to equation (0; however the electronic 
part, independent of a:, is not the clamped-nuclei Hamiltonian. 
Instead, the exact electronic Hamiltonian can be expressed as 
a direct integral of clamped-nuclei Hamiltonians and necessar- 
ily has a purely continuous spectrum of energy levels; there are 
no potential energy surfaces. This continuum has nothing to do 
with the molecular centre-of-mass, by construction. The paper 
concludes (® with a discussion of our finding. 

2. The Born-Oppenheimer approximation 

The original Born and Oppenheimer approximation is 
summarized in the famous book by Born and Huang, and the 
later Born adiabatic method [6] is given in an appendix to that 
book [7]. Born and Huang use the same notation for both for- 
mulations and it is convenient to follow initially their presen- 
tation; the following is a short account focusing on the main 
ideas. They work in a position representation and for simplicity 
suppress all individual particle labels. Let us consider a system 
of electrons and nuclei and denote the properties of the former 
by lower-case letters (mass m, coordinates x, momenta p) and 
of the latter by capital letters (mass M, coordinates X, momenta 
P). The kinetic energy of the nuclei is the operator 

n 2 ( d 2 
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2M 



2M\dX 2 



and that of the electrons 

2m 



P 2 = 



l—l 2m \ dx 2 



(1) 



(2) 



The total Coulomb energy of the electrons will be represented 
by U(x, X). We further introduce the abbreviation 



+ U = H \x,—,X 
ox 



(3) 



The full Hamiltonian for the molecule is then 



H = T e + U +T N = H„ + T N 



(4) 



The fundamental idea of Born and Oppenheimer is that the 
low-lying excitation spectrum of a typical molecule can be cal- 
culated by regarding the nuclear kinetic energy TV as a small 
perturbation of the Hamiltonian H . The physical basis for the 
idea is the large disparity between the mass of the electron and 
all nuclear masses. The expansion parameter must clearly be 
some power of m/M , where M a can be taken as any one of 
the nuclear masses or their mean. They found that the correct 
choice is 

m 



and therefore 



k 4 H x 



d_ 

dX 



it 

2m 



dX 2 



(5) 



Thus the total Hamiltonian may be put in the form 

H = H + k 4 H x (6) 

with Schrodinger equation 

(H - E)i//(x,X) = 0. (7) 

In the original paper Born and Oppenheimer say at this point 
in their argument that [ 3 ] : 

If one sets k — ... one obtains a differential equa- 
tion in the x alone, the X appearing as parameters: 



ij/ = 0. 



This represents the electronic motion for stationary 
nuclei. 

and it is perhaps to this statement that the idea of an elec- 
tronic Hamiltonian with fixed nuclei as arising by letting the 
nuclear masses increase without limit, can be traced. In modern 
parlance H is customarily referred to as the "clamped-nuclei 
Hamiltonian". 

Consider the unperturbed electronic Hamiltonian H„(x,Xf) 
at a fixed nuclear configuration Xf that corresponds to some 
molecular structure. The Schrodinger equation for H is 



(H (x,X f ) - E°(X f ) m y(x,X f ) m = 0. 



(8) 



This Hamiltonian's natural domain, D„, is the set of square in- 
tegrable electronic wavefunctions [<p m } with square integrable 
first and second derivatives; T> is independent of Xf. We may 
suppose the {<p m } are orthonormalized independently of Xf 



I 



dxtp(x,Xf)* n <p(x,Xf)„ 



In the absence of degeneracies ("curve-crossing") they may be 
chosen to be real; otherwise there is a phase factor to be consid- 
ered. For every Xf, H a is self-adjoint on the electronic Hilbert 
space 'H(Xf), and therefore the set of states {<p(x,Xf) m } form a 
complete set for the electronic Hilbert space indexed by Xf. 

The clamped-nuclei Hamiltonian can be analysed with the 
HVZ theorem which shows that it has both discrete and contin- 
uous parts to its spectrum ifs - 121. 



a(X f )= a{H {x,X f )) = 



E°(X f ) , . . . E°(X f )Jj (J [a(X/), oo) 

(9) 

where the {E°(Xf)i c } are isolated eigenvalues of finite multiplic- 
ities. A(Xj ) is the bottom of the essential spectrum marking the 
lowest continuum threshold. In the case of a diatomic molecule 
the electronic eigenvalues depend only on the internuclear sep- 
aration r, and have the form of the familiar potential curves 
shown in FigQ] For the general polyatomic molecule, the dis- 
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Figure 1: The spectrum cr(r) for a diatomic molecule | 8]. 

crete eigenvalues are molecular potential energy surfaces. 

Born and Oppenheimer used the set {ip m } to calculate ap- 
proximate eigenvalues of the full molecular Hamiltonian H on 
the assumption that the nuclear motion is confined to a small 
vicinity of a special (equilibrium) configuration X®. Their great 
success was in establishing that the energy levels of the low- 
lying states typical of small polyatomic molecules obtained 
from molecular spectroscopy could be written as an expansion 
in powers of k 2 
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V< 0) 



+ K' 



P(2) 



K*E 



(4) 
nvJ 



(10) 



where V„ is the minimum value of the electronic energy 

(2) 

which characterises the molecule at rest, E m ; is the energy of 
the nuclear vibrations, and E^J, contains the rotational energy 
12]. The corresponding approximate wavefunctions are simple 
products of an electronic function <p m and a nuclear wavefunc- 
tion; this is known as the adiabatic approximation. In the orig- 
inal perturbation formulation the simple product form is valid 
through k 4 , but not for higher order terms. 

About 25 years later Bom observed that the results of molec- 
ular spectroscopy suggest that the adiabatic approximation has 



a wider application than predicted by the original theory, and he 
proposed an alternative formulation . It is assumed that the 
functions E°(X) m and <p(x, X) m arising from equation © which 
represent the energy and wavefunction of the electrons in the 
state m for a fixed nuclear configuration X, are known. Born 
proposed to solve the wave equation (I7J by an expansion 



(ii) 



with coefficients {0(X) m } that play the role of nuclear wave- 
functions. Substituting this expansion into the full Schrodinger 
equation Q, multiplying the result by <p(x, X)* and integrating 
over the electronic coordinates x leads to a system of coupled 
equations for the nuclear functions {O}, 



(T N + E°(X)„ 



E)Q>(X) n + YjC(X,P)„ 



.<B(X)„, = (12) 



where the coupling coefficients \C(X, P) ml >) have a well-known 
form which we need not record here |]7|. In this formulation the 
adiabatic approximation consists of retaining only the diagonal 
terms in the coupling matrix C(X, P), for then 



\f/{x,X) 



<p(x,X) n ®(X) n . 



(13) 



An obvious defect in this presentation, recognized by Born 
and Huang [7], is that the kinetic energy of the overall center- 
of-mass is retained in the Schrodinger equation (IT); this is easily 
corrected, either by the explicit separation of the center-of-mass 
kinetic energy operator (see $3), or implicitly, as in the compu- 
tational scheme proposed by Handy and co-workers fl 1 3l — 1 1 3h - 
For many years now these equations have been regarded in the 
theoretical molecular spectroscopy/quantum chemistry litera- 
ture as defining the "Bom-Oppenheimer approximation", the 
original perturbation method being relegated to the status of his- 
torical curiosity. Commonly they are said to provide an exact 
(in principle) solution II 16142011 for the stationary states of the 
molecular Schrodinger equation (U_), it being recognized that in 
practice drastic truncation of the infinite set of coupled equa- 
tions (fT2l is required. In the next section we make a critical 
evaluation of this conventional account. 



3. The molecular Schrodinger equation 

We now start again and reconsider the Hamiltonian for a col- 
lection of electrons and nuclei, assuming their interactions are 
restricted to the usual Coulombic form. The key idea in §|2] 
is the decomposition of the molecular Hamiltonian <j4j into a 
part containing all contributions of the nuclear momenta, and 
a remainder. This must be done in conjunction with a proper 
treatment of the center-of-mass motion. These two ideas guide 
the following discussion. 

Let the position variables for the particles in a laboratory 
fixed frame be designated as fx, }. When it is necessary to distin- 
guish between electrons and nuclei, the position variables may 
be split up into two sets, one set consisting of N coordinates, 
x e , describing the electrons with charge —e and mass m, and 



3 



the other set of A coordinates, x™, describing the nuclei with 
charges +Z g e and masses m g , g — I,.,. A; then n — N + A. We 
define corresponding canonically conjugate momentum vari- 
ables {p,} for the electrons and nuclei. With this notation, and 
after the usual canonical quantization, the Hamiltonian operator 
for a system of N electrons and A atomic nuclei with Coulombic 
interactions may be written as 



H 



A 



+ y £L _ _f_ y ti 



(14) 



It is easily seen that the Coulomb interaction is translation 
invariant. Thus the total momentum operator, which in the 
present notation is 

n 
i 

commutes with H. Physically the center-of-mass of the whole 
system, with position operator 

j n n 

R = ^:Z m ' x " M T = Yu m 

i i 

behaves like a free particle. It is desirable then to introduce R 
and its conjugate Pr, together with appropriate internal coordi- 
nates, into H to make explicit the separation of the center-of- 
mass and the internal dynamics. Formally H may be written as 
a direct integral [21] 



H 



Jr 3 



where [22] 



H(P) 



H(P) dP 



+ H' 



2M T 



(15) 



(16) 



is the Hamiltonian at fixed total momentum P. This representa- 
tion shows directly that H has purely continuous spectrum. The 
internal Hamiltonian H' is independent of the center-of-mass 
variables and acts on L 2 (R 3( " _1) ); it is explicitly translation in- 
variant. 

There are infinitely many possible choices of internal coor- 
dinates that are unitarily equivalent, so that the form of H' is 
not determined uniquely, but whatever coordinates are chosen 
the essential point is that H' is the same operator specified by 
the decomposition (I161 l. A simple procedure to make the inter- 
nal Hamiltonian explicit is to refer the particle coordinates to a 
point moving with the system, for example, the center-of-mass 
itself, the center-of-nuclear-mass or one of the moving particles 
ill. 

As a result of the transformation to internal variables the ki- 
netic energy operators are no longer diagonal in the particle 
indices and certain choices of the moving point, such as the 
choice of a single nucleus, result in an operator in which the 
nuclear and electronic indices are mixed. However the choice 
of the center-of-nuclear-mass as the point of origin avoids this 



mixing; the explicit equations for this choice were given in Sut- 
cliffe and Woolley 12411 . There are A - 1 translationally invariant 
coordinates t" expressed entirely in terms of the original x^j that 
may be associated with the nuclei, and there are N translation- 
ally invariant coordinates for the electrons t? which are simply 
the original electronic coordinates referred to the center-of- 
nuclear-mass. There are corresponding canonically conjugate 
internal momentum operators. The classical total kinetic en- 
ergy still separates in the form 

To = Tn + T e = T C M + Tjv h + J e i 

and the same is true after quantization. 1cm ls tne kinetic en- 
ergy for the center-of-mass and, for example, T e / only involves 
electronic variables 



1=1 u=i 



f) • <«5> 



(17) 



with 



M 



111 
m M 



The nuclear kinetic energy term is like the second term in 
( TPTI i expressed in terms of internal nuclear momentum variables 
{II"} with a mass factor in the denominator composed from the 
nuclear masses. 

After this unitary transformation the original Coulomb 
Hamiltonian operator for the molecule can be rewritten in the 
form 

H(x e ,p e ,x'\p n ) = 1 CM + J Nu + H elec (18) 

where H e ^ ec is composed of ( fT77T > together with all the Coulomb 
interaction operators expressed in terms of the jt e ,t n ) position 
operators. 

The spectrum of the Coulomb operator H' formed by drop- 
ping J cm from ( TT81 was considered by Kato ll26Tl who showed 
in Lemma 4 of his paper that for a Coulomb potential V and 
for any function / in the domain Do of the full kinetic energy 
operator To, the domain, Oy, of the internal Hamiltonian H' 
contains Do and there are two constants a, b such that 

l|V/||<a||T /||+ 6H/II 

where a can be taken as small as is liked. This result is of- 
ten summarised by saying that the Coulomb potential is small 
compared to the kinetic energy. Given this result he proved in 
Lemma 5 (the Kato-Rellich theorem) that the usual operator is 
indeed, for all practical purposes, self-adjoint and so is guar- 
anteed a complete set of eigenfunctions, and is bounded from 
below. 

In the present context the important point to note is that the 
Coulomb term is small only in comparison with the kinetic 
energy term involving the same set of variables. So the ab- 
sence of one or more kinetic energy terms from the Hamiltonian 
may mean that the Coulomb potential term cannot be treated as 
small. Such a defective Hamiltonian will no longer be self- 
adjoint in the way demonstrated by Kato. 
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With our choice of coordinates the translationally invariant 
Hamiltonian may be written as the sum of the last two terms in 
( fl~8b and the electronic Hamiltonian now becomes 



elec 



& JL 



2m 



!=1 

2 A N 



Ane< 



yy^_ 



where it is understood that the x" are to be realised by a suit- 
able linear combination of the tf. The electronic Hamiltonian is 
properly translationally invariant and, were the nuclear masses 
to increase without limit, the second tem in $1% would vanish 
and the electronic Hamiltonian would have the same form as 
that usually invoked in molecular electronic stucture calcula- 
tions. If the nuclear positions were chosen directly as a trans- 
lationally invariant set, it would be those values which would 
appear in the place of the nuclear variables. The nuclear part 
involves only kinetic energy operators and has the form: 



I Nu 



Cm "J 



(20) 



with the inverse mass matrix fi given in terms of the coefficients 
that define the translationally invariant coordinates and the nu- 
clear masses. 



H 



elec 



has the same invariance under the rotation-reflection 



group 0(3) as does the full translationally invariant Hamilto- 
nian H' and, in the case of the molecule containing some iden- 
tical nuclei, it has somewhat extended invariance under nuclear 
permutations, since the nuclear masses appear only in symmet- 
rical sums. In a position representation the Schrodinger equa- 
tion for H e ^ ec may be written as 



(21) 



where m is used to denote a set of quantum numbers (JMprk): 
J and M for the angular momentum state: p specifying the par- 
ity of the state: r specifying the permutationally allowed irreps 
within the group(s) of identical particles and k to specify a par- 
ticular energy value [24]. 

Before discussing (l2"TT i for the general case, let us consider a 
specific example for which high-quality computational results 
exist. It might be hoped, in the light of the claim in the origi- 
nal paper by Born and Oppenheimer quoted in $2] that the so- 
lutions of (|2"T1) would actually be those that would have been 
obtained from an exactly solved problem by letting the nuclear 
masses increase without limit. Although there are no exactly 
solved molecular problems, the work of Frolov [25] provides 
extremely accurate numerical solutions for a problem with two 
nuclei and a single electron. Frolov investigated what happens 
when the masses of one and then two of the nuclei increase 
without limit in his calculations. To appreciate his results, con- 
sider a system with two nuclei; the natural nuclear coordinate is 



the internuclear distance which will be denoted here simply as 
t. When needed to express the electron-nuclei attraction terms, 
x n is simply of the form or,t where cr, is a signed ratio of the 
nuclear mass to the total nuclear mass; in the case of a homonu- 
clear system a-, = ±A. 

The di-nuclear electronic Hamiltonian after the elimination 
of the center-of-mass contribution is 



2m 

1=1 

9 N 

- — Y 

4/reo j-( 

2 N 



2{m\ + m 2 ) 
Z 2 



Y v(t?) • v(tp 



t e + ant 



e 1 y ' 1 



It) + a 2 t\ ) 
ZiZ 2 



R 



R = 



while the nuclear kinetic energy part is: 

Tv„(t) = -^(- + -)v 2 (t). 

L \m\ mi) 



2/j 



V 2 (t). 



(22) 



(23) 



The full internal motion Hamiltonian for the three-particle sys- 
tem is then 

H'(t e ,t) = H' e (t e ) + Tv„(t). (24) 

It is seen from (l23l . that if only one nuclear mass increases 
without limit then the kinetic energy term in the nuclear variable 
remains in the full problem and so the Hamiltonian (l24l remains 

self-adjoint in the Kato sense. Frolov's calculations showed 
that when one mass increased without limit (the atomic case), 
any discrete spectrum persisted but when two masses were al- 
lowed to increase without limit (the molecular case), the Hamil- 
tonian ceased to be well-defined and this failure led to what 
he called adiabatic divergence in attempts to compute discrete 
eigenstates of d24l i. This divergence is discussed in some math- 
ematical detail in the Appendix to Frolov(25J. It does not arise 
from the choice of a translationally invariant form for the elec- 
tronic Hamiltonian; rather it is due to the lack of any kinetic 
energy term to dominate the Coulomb potential. Thus it really 



1 elec 



to see whether 



is essential to characterize the spectrum of H 
the traditional approach can be validated. 

As before, the nuclear kinetic energy operator is proportional 
to k 4 , so after dropping the uninteresting center-of-mass kinetic 
energy term, (fT¥t is seen to be of the same form as ©. There is 
however a fundamental difference between © and ([T8l > which 
may be seen as follows; with the center-of-nuclear-mass chosen 
as the electronic origin H e ^ ec is independent of the nuclear mo- 
mentum operators and so it commutes with the nuclear position 
operators 

[H elec ,t n ] = 0. 

They may therefore be simultaneously diagonalized and we use 
this property to characterize the Hilbert space H for H e ' ec . 
Let b be some eigenvalue of the t n corresponding to choices 
[x ? = a g ,g = I,... A) in the laboratory-fixed frame; then 
the {a g } describe a classical nuclear geometry. The set, X, of 

all b is R 3(/i_I) . We denote the Hamiltonian H elec evaluated 
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at the nuclear position eigenvalue b as H(b, t e ) = H for 
short; this H c is very like the usual clamped-nuclei Hamilto- 
nian but it is explicitly translationally invariant, and has an extra 
term, the second term in (TTTl i. (or (l22t). which is often called 
the Hughes-Eckart term, or the mass polarization term. The 
Schrodinger equation for H is of the same form as (H), with 
eigenvalues E°(b)i l and corresponding eigenfunctions y>(t e ,b)*, 
and with spectrum <x(b) analogous to 

In general the symmetries of H(b, t e )„ considered as a func- 
tion of the electronic coordinates, t e , are much lower than those 
of H e ^ ec , since they are determined by the (point group) trans- 
formations that leave the geometrical figure defined by the {a,,} 
invariant. In a chiral structure there is no symmetry operation 
other than the identity so even space-inversion symmetry (par- 
ity) is lost. 

At any choice of b the eigenvalues of H„ will depend only 
upon the shape of the geometrical figure formed by the {a s } 
and not at all upon its orientation. For other than diatomic 
molecules, it is technically possible to produce a new coordi- 
nate system from the b with angular variables describing the 
orientation of the figure and 3A-6 internal variables describ- 
ing its geometry. However this cannot be done without angular 
momentum terms arising in which electronic and nuclear mo- 
mentum operators are coupled. Thus the kinetic energy can no 
longer be expressed as a sum of separate electronic and nuclear 
parts. There are also many possible choices of angular and in- 
ternal variables which span different (though often overlapping) 
domains, so it is not possible to produce a single account of the 
separation of electronic and nuclear motion, applicable to all 
choices. 

As in §|2] H e is self-adjoint on an electronic Hilbert space 
"H(b), so we have a family of Hilbert spaces { r H(b)} which are 
parameterized by the nuclear position vectors b e X that are the 
"eigenspaces" of the family of self-adjoint operators H ; from 
them we can construct a big Hilbert space as a direct integral 
over all the b values 

"H = H(b)db (25) 
Jx 

and this is the Hilbert space for H elec in COD. The internal 
molecular Hamiltonian H' in (TToT l and the clamped-nuclei like 
operator H„ just defined can be shown to be self-adjoint (on 
their respective Hilbert spaces) by reference to the Kato-Rellich 
theorem ||9|] because in both cases there are kinetic energy op- 
erators that dominate the (singular) Coulomb interaction; they 
therefore have a complete set of eigenfunctions. This argu- 
ment cannot be made for H e ^ ec because it contains nuclear 
position operators in some Coulombic terms but there are no 
corresponding nuclear kinetic energy terms to dominate those 
Coulomb potentials. It therefore cannot be thought of as being 
self-adjoint on the full Hilbert space of electronic and nuclear 
variables. Thus even though H e ' ec is evidently Hermitian sym- 
metric, this is not sufficient to guarantee a complete set of eigen- 
functions because the operator is unbounded, and an expansion 
analogous to (ITTb seems quite problematic. 

However that may be, equation (1231 1 leads directly to a fun- 



damental result; since H commutes with all the {t n }, it has 
the direct integral decomposition 

H elec _ f H (b,t e ) db. (26) 
Jx 

This result implies at once that the spectrum of H e ' ec is purely 
continuous 

o- = <r(H elec ) = y o-(b) = [Vo, oo) 

b 

where Vo is the minimum value of E(b)o\ in the diatomic 
molecule case this is the minimum value of £b(r) (see Fig. [TJ. 
I_|elec nas no norma ij za ble eigenvectors. We conclude there- 
fore that the decomposition of the molecular Hamiltonian into 
a nuclear kinetic energy operator contribution, proportional to 
k 4 , and a remainder does not yield molecular potential energy 
surfaces. 

4. Discussion 

The rigorous mathematical analysis of the original pertur- 
bation approach proposed by Born and Oppenheimer (]2| for 
a molecular Hamiltonian with Coulombic interactions was ini- 
tiated by Combes and co-workers ||8l \2Jl |28|] with results for 
the diatomic molecule. A perturbation expansion in powers 
of k leads to a singular perturbation problem because k is 
a coefficient of differential operators of the highest order in 
the problem; the resulting series expansion of the energy is 
an asymptotic series, closely related to the WKB approxima- 
tion. Some properties of the operator H e ^ ec , (equation [26), 
seem to have been first discussed in this work. Since the ini- 
tial work of Combes, a considerable amount of mathematical 
work has been published using both time-independent and time- 
dependent techniques with developments for the polyatomic 
case; for a recent review of rigorous results about the separa- 
tion of electronic and nuclear motions see Hagedorn and Joye 
lE9ll which covers the literature to 2006. 

If it can be assumed that a) the electronic wavefunction van- 
ishes strongly outside a region close to a particular nuclear ge- 
ometry and b) that the electronic energy at the given geome- 
try is an isolated minimum, then it is possible to present a rig- 
orous account of the separation of electronic and nuclear mo- 
tion which corresponds in some measure to the original Born- 
Oppenheimer treatment; such an account is provided for a poly- 
atomic molecule by Klein and co-workers II3CH1 . Because of the 
continuous spectrum of the electronic Hamiltonian H e ^ ec , it is 
not possible to use regular perturbation theory in the analysis; 
instead asymptotic expansion theory is used so that the result 
has essentially the character of a WKB approximation. Nei- 
ther is it possible to consider rotational symmetry for the rea- 
sons outlined in §|3] However it is possible to impose inversion 
symmetry and the nature of the results derived under this re- 
quirement depend on the shape of the geometrical figure at the 
minimum energy configuration. If the geometry at the mini- 
mum is either linear or planar then inversion can be dealt with 
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in terms of a single minimum in the electronic energy. If the ge- 
ometry at the minimum is other than these two forms, inversion 
produces a second potential minimum and the problem must be 
dealt with as a two-minimum problem; then extra consideration 
is necessary to establish whether the two wells have negligi- 
ble interaction so that only one of the wells need be considered 
for the nuclear motion. The nuclei are treated as distinguish- 
able particles that can be numbered uniquely. The symmetry 
requirements on the total wavefunction that would arise from 
the invariance of the Hamiltonian operator under the permuta- 
tion of identical nuclei are not considered. 

It is similarly possible to consider such phenomena as 
Landau-Zener crossing by using a time-dependent approach to 
the problem and looking at the relations between the electronic 
and nuclear parts of a wave packet 03111 . This is essentially a 
use of standard coherent state theory where again the nuclei are 
treated as distinguishable particles and the method is that of 
asymptotic expansion. 

To summarize: what the rigorous work so far presented has 
done is to show that if an electronic wavefunction has certain 
local properties and if the nuclei can be treated as distinguish- 
able particles then the the eigenfunctions and eigenvalues of the 
molecular (Coulomb) Hamiltonian can be obtained in a WKB 
expansion in terms of k to arbitrary order. 

In this paper we have attempted to discuss the Born approach 
to the molecular theory in terms first set out by Combes n27l . 
The essential point is that the decomposition of the molecular 
Hamiltonian (with center-of-mass contribution removed) into 
the nuclear kinetic energy, proportional to k 4 and a remainder, is 
specified by equation (fist , not by (0, or in other words, equa- 
tion (|6|) cannot be written with an — sign. Allowing the nuclear 
masses to increase without limit in H e ^ ec does not produce an 
operator with a discrete spectrum since this would just cause 
the mass polarisation term to vanish and the effective electronic 
mass to become the rest mass. It is thus not possible to re- 
duce the molecular Schrodinger equation to a system of cou- 
pled differential equations of classical type for nuclei moving 
on potential energy surfaces, as suggested by Born, without a 
further approximation of an essentially empirical character. An 
extra choice of fixed nuclear positions must be made to give any 
discrete spectrum and normalizable L 2 eigenfunctions. In our 
view this choice, that is, the introduction of the clamped-nuclei 
Hamiltonian, by hand, into the molecular theory as in ^2]is the 
essence of the "Born-Oppenheimer approximation" 



H 



elec 



I H(b,f) db -> H(b,t e )„ 

Jx 

Born-Oppenheimer approximation. (27) 



If the molecular Hamiltonian H were classical, the removal of 
the nuclear kinetic energy terms would indeed leave a Hamil- 
tonian representing the electronic motion for stationary nuclei, 
as claimed by Born and Oppenheimer [4 |3J. As we have seen, 
quantization of H changes the situation drastically, so an im- 
plicit appeal to the classical limit for the nuclei is required. The 
argument is a subtle one, for subsequently, once the classical 
energy surface has emerged, the nuclei are treated as quantum 



particles (though indistinguishability is rarely carried through); 
this can be seen from the complexity of the mathematical ac- 
count given by Klein and co-workers l30ll . 

We have long argued that the solutions of the time- 
independent Schrodinger equation for the molecular Hamilto- 
nian are of limited interest for chemistry, being really only rel- 
evant to a quantum mechanical account of the physical prop- 
erties (mainly spectroscopic) of atoms and diatomic molecules 
in the gas-phase. Towards the end of his life, P.-O. Lowdin 
made an extended study of a quantum mechanical definition of 
a molecule; in one of his late papers 113211 he lamented 



The Coulombic Hamiltonian H' does not provide 
much obvious information or guidance, since there 
is [sic] no specific assignments of the electrons oc- 
curring in the systems to the atomic nuclei involved - 
hence there are no atoms, isomers, conformations etc. 
In particular one sees no molecular symmetry, and 
one may even wonder where it comes from. Still it 
is evident that all this information must be contained 
somehow in the Coulombic Hamiltonian. 

In our view it is not at all evident that the Coulombic Hamil- 
tonian on its own will give rise to the chemically interesting 
features Lowdin required of it, nor will they be approachable 
by regular perturbation theory (supposedly convergent) starting 
from its eigenstates. Fundamental modifications of the quan- 
tum theory of the Coulomb Hamiltonian for a generic molecule 
have to be made for a chemically significant account of dipole 
moments, functional groups and isomerism, optical activity and 
so on I33I437I1 . In other words one should not expect useful con- 
tact between the quantum theory of an isolated molecule (which 
is what the eigenstates of the Coulombic Hamiltonian refer to) 
and a quantum account of individual molecules, as met in or- 
dinary chemical situations where persistent interactions (due 
to the quantized electromagnetic field, other molecules in bulk 
media) and finite temperatures are the norm [38, 39]. 

The usual practice in computational chemistry in which 
clamped-nuclei electronic energy calculations are used to define 
a potential energy surface upon which a quantum mechanical 
nuclear motion problem is solved, is a practice that is well de- 
fined mathematically. Solutions obtained to the equations for- 
mulated in this context do not have the symmetry of the full 
internal motion Hamiltonian and their relationship to solutions 
of that Hamiltonian is, as has been seen, uncertain. 

That said, the conventional account, treating formally iden- 
tical nuclei as identifiable particles when it seems chemically 
prudent to do so, has enabled a coherent and progressive ac- 
count of much chemical experience to be provided. But it is not 
derived by continuous approximations from the eigensolutions 
of the Schrodinger equation for the molecular Coulomb Hamil- 
tonian, requiring as it does an essential empirical input. The 
tremendous success of the usual practice might perhaps be best 
regarded as a tribute to the insight and ingenuity of the practi- 
tioners for inventing an effective variant of quantum theory for 
Chemistry. 
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